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SUMMARY 

In 3une 1978 and July 1979, the Air Resources Branch carried out intensive 
field programs in the Sudbury area to study the dispersion of the Inco 381 metre 
stack plume in the convective boundary layer. The studies were designed to 
provide experimental data for the verification of a short term, short range 
fumigation model (it) for the estimation of ground level concentrations of S0 2 due 
to this stack's emissions. 

The studies provided a set of data suitable for checking the validity of the 
model predictions. They also provided a set of intense meteorological 
measurements from which characteristics of the convective boundary layer in the 
Sudbury area were determined. 

It was found that the convective boundary layer exhibits a typical pattern 
under relatively warm sunny conditions. It grows rapidly from roughly 0900 until 
roughly 1300 hours local time after which it remains relatively constant until 
shortly before sunset. 

The height of the convective layer is generally demarked by a strong 
capping inversion. Within the layer (in a barotropic CBL) the wind speed and 
potential temperature is approximately constant with height. The relatively long- 
lived convective turbulence elements within the layer bring about rapid vertical 
mixing which often causes portions of the plume to reach ground level very close to 
the stack. 
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CONCLUSIONS 

Field studies were carried out in June 1978 and August 1979 into dispersion 
of the INCO 381 m stack plume in the convective boundary layer (CBL). The 
studies provided data for the validation of a computer model for estimating 
ground level SO ? concentrations resulting from the INCO 381 M stack emissions 
under convective conditions (4). The data are reported in a format directly 
applicable to this purpose. 

The model development and validation are published elsewhere (3,4). A good 
comparison was obtained between the model predictions and the field study data. 

In addition to providing data for model validation, the field studies also 
provided an opportunity to study the characteristics of the convective boundary 
layer in the Sudbury area. A number of conclusions resulted from this work. They 
are: 

• The CBL develops from roughly 0900 to 1 300 hours local 
time after which it remains quasi-steady-state for most 
of the afternoon . In the evening, the dissolution of the 
CBL is very rapid and occurs just before sunset. 

• The mixed layer is generally characterized by a capping 
inversion, a constant wind speed and a constant potential 
temperature lapse rate (equivalent to the dry adiabatic 
lapse rate). 

• Plume dispersion within the CBL is characterized by 
behaviour known as looping. On warm sunny days the 
INCO Supserstack plume begins looping typically 
between 1000 and 1100 hours and stops roughly one to 
two hours before sunset. The initial impingment of 
individual plume loops usually occurs at a distance of 3 
to 10 km from the stack and results in intermittent, 
high, ground level concentrations. Farther downwind, 
the plume becomes well-mixed in the vertical. 
Concentrations are not as high, and less intermittent. 



1. INTRODUCTION 

In Dune 1978 and August 1979, the Air Resources Branch of the Ministry to 
the Environment carried out two intensive field studies in the Sudbury area. These 
studies were designed to investigate dispersion of the INCO 381 m stack piume in 
the convective boundary layer. Both studies were carried out under the auspices of 
the Sudbury Environmental Study. 

The objectives of these intensive field studies are related to the very nature 
of dispersion of elevated pollutant releases in the convective boundary layer. A 
brief description of the latter is therefore given here based on other publications 
(1-3). The reader is referred to these publications for more detailed descriptions. 

The boundary layer over land can be thought of as a three layer structure. 
Each of the layers is characterized by a separate set of parameters related to the 
turbulence in it. Beginning at the surface, the layers are: 

1. The Surface Layer: 

In the surface layer, wind shear is dominant in the production of turbulence. 
The layer is confined to a height range Z< I l| where Z is height above ground and 
L is defined as the Obukhov length (=-u 3 */k(g/T)H ; with u* being the friction 
velocity, T the average temperature in the layer, H the kinematic surface heat 
flux and k the von Karman constant). Within the layer the controlling parameters 
are Z, T , H and g/T where T is the surface shear stress. These parameters can be 
combined into a scaling velocity and temperature respectively , i.e., 

y* = (T / P )K (1) 

T* = -Ho/u* (2) 



Dimensionless groups formed with U* and T* become universal functions of 
Z/L within the layer. 

2. The Free Convection Layer: 

In this layer, shear stress is no longer dominant but height, Z, is still the 
significant length scale. It is suggested (1) that the upper limit of this sublayer is 
0.1 Z\ where Zj is the height of the mixed layer (discussed next). In the free 
convection layer the governing parameters are Z\, H , and g/T and the scaling 
velocity and scaling temperature are: 



Uf= <H Z(g/T)) i/3 (3) 

T f = H /u f (*) 



3. The Mixed Layer: 

The mixed layer is characterized by turbulence insensitive to height Z and 
shear stress x . The height of the mixed layer is defined as the height of the 
lowest inversion base Z\ capping the layer and this height in the controlling length 
scale. The scaling velocity and scaling temperature are: 

w* = (H Zi (g/T)) 1/3 (5) 

9* = H /w* (6) 



Within the mixed layer, dimensionless groups formed with w* and 0* should 
be functions of Z/Z[ only. 

The term "mixed layer" is generally used synonymously with the term 
"convective boundary layer (CBL)" which is defined as the portion of the 
atmosphere most directly affected by solar heating of the earth's surface. The 
heating results in the buoyant production of turbulence in the form of long-lived 



updrafts and downdrafts. Rapid vertical mixing takes place because of these 
motions and this mixing extends to the top of the layer - generally one to two 
kilometers. 

Under conditions of very strong mixing, the CBL is generally barotropic, 
i.e., surfaces of constant pressure are also surfaces of constant density. In these 
cases, wind speed does not change with height and the vertical profiles of both 
wind speed and potential temperature are constant with height to the top of the 
layer. 

Smoke plumes emitted from elevated sources within the CBL exhibit 
behaviour described rather graphically as 'looping'. Looping results from portions of 
the plumes being caught in alternating, large-scale, turbulent updrafts and 
downdrafts immediately after release from the stack. Those portions of plumes 
emitted into downdrafts are often carried to ground level very close to the stack 
causing high ground level concentrations. Those portions of plumes emitted into 
updrafts are carried high into the mixed layer resulting in zero ground level 
concentrations near the stack. Further downwind these portions do reach ground 
level causing non-zero ground level concentrations. 

Downwind of the stack, the strong vertical motions tend to mix the 
pollutants rather rapidly throughout the entire mixed layer and pollutant 
concentrations are relatively homogeneous in the vertical. 

This report describes two intensive field studies carried out in the Sudbury 
area to investigate the high ground level concentrations caused by the INCO 381 m 
stack plume under convective conditions. This behaviour has been observed to 
occur from mid-April to late October on sunny days. 
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The Air Resources Branch under the auspices of the Sudbury Environmental 
Study decided to investigate this phenomenon and to develop a computer model to 
calculate impingement points and ground level concentrations when the plume is 
emitted into the CBL. Ultimately the model will be used to assess the impact of 
local emissions on air quality and dry deposition in the Sudbury area. 

The computer model was developed by the Air Quality and Meteorology 
Section of the Air Resources Branch. To aid in the development and verification of 
this model, the Technology Development and Appraisal Section (now the 
Atmospheric Research and Speical Programs Section) planned and carried out both 
field studies. The experimental portions of the study were designed to 
simultaneously measure: (1) all input variables to the model (2) the location and 
magnitude of ground level concentrations of sulfur dioxide emitted from the INCO 
381 m stack and (3) characteristics of the CBL in Sudbury. 

This report describes the instrumentation, field programs and results of both 
studies. The reader is referred to other publications (3,4) for a detailed 

description of the model and its comparison with the field study data. 



2. EXPERIMENTAL PROGRAM 

The experimental programs for both field studies were designed to provide 
data for the verification of the computer model developed by the Air Quality and 
Meteorology Section and to study the characteristics of the convective boundary 
layer in the Sudbury area. In order to do this, three types of measurements were 
undertaken simultaneously, namely: 



1. Ground level concentration measurements of SOj taken for 
comparison with predictions of the model. 

2. Meteorological measurements of total incoming solar radiation, S 
(W/m 2 ); mixed layer wind speed, U (m/s); ambient temperature at 
stack height, T (°C); and mixed layer height, Z. (m). Data from 
these and surface heat flux measurements were also used to 
characterize the structure of the CBL. 

3. In-stack measurements of S0 2 emission rate, Q g (g/s); stack gas 

temperature, T (°C); and stack gas velocity, V (m/s). These 
s ** 

measurements were carried out by INCO as part of their Emission 
Reduction Program and data were provided to the Ministry of the 
Environment after each study. These measurements are to be used 
as input to the model. 

The instrumentation and methods used for these measurements are 
discussed in detail below. A brief summary is also given in Tables 1(a) and Kb). 

2.1 GROUND LEVEL CONCENTRATION MEASUREMENTS 

The instrumentation for the ground level concentration measurements was 
similar for both field studies. It consisted of two measurement platforms- each 
outfitted with a fast-response S0 2 continuous analyser - the Sign -X, a Hewlett- 
Packard Strip Chart Recorder and a Topaz 1 Kw inverter. In the June study, one 
platform was mounted in a van and the other in a twin-engine Piper Aztec aircraft. 
In the August study, both platforms were mounted in station wagons and no aircraft 
was used. In all cases, the inlet to the Sign - X was mounted outside the region of 
aerodynamic distortion about the transport unit. 
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The vehicles used in the August 1979 study carried additional support 
instrumentation for assisting in the concentration measurements. Both carried 
ground-recovery intervalometers which maintained the speed of the chart 
recorders at a speed dependent on that of the vehicle. One of the vehicles 
(instrumented by and rented from Moniteq Ltd. of Toronto), also carried a MARSS 
100 UV Absorption Spectrometer. This instrument was a prototype manufactured 
by Moniteq Ltd. and was taken to Sudbury for field testing. It was used to locate 
the overhead burden of S0 2 which helped define the boundaries of the plume. 

A Barringer Research Ltd Correlation Spectrometer (COSPEC) was also 
provided by Moniteq and was used in the other vehicle for the same purpose. Use 
of the latter type of spectrometer for plume tracking has been well documented 
(5,6) and both instruments proved very useful in this study. 

The experimental program was designed to provide simultaneous ground 
level SO- concentration measurements at two different distances downwind of the 
stack. One Sign - X was therefore deployed in the near field (5 to 15 Km from the 
stack) and the other in the far field (15 to 40 Km). The aircraft was always used in 
the far-field because air transport regulations prevented low level flying over areas 
near the stack. It should be pointed out that it was not always possible to obtain 
simultaneous near-field and far-field measurements because of logistical problems. 
However, the data from one instrument only was also satisfactory for the model 
verification. 

The field measurement procedures were relatively simple. Each day, at the 
onset of convective activity, the plume direction was noted. A topographical map 
was used to locate roadways running roughly perpendicular to the axis of the 
plume. The vehicles were deployed to these roads where they 



carried out repeated cross-wind traverses roughly perpendicular to the plume axis. 
Naturally, the traverses followed existing roadways so that the directions of 
individual portions of the traverses were variable. Bench marks were located along 
each traverse path and marked on the chart record of the Sign - X ouput. This 
allowed easy ground path recovery for later data analysis. 

Five to twenty cross-wind traverses were made at each location within the 
duration of an hour. The importance of this number and time of traverses is 
discused in another publication (7). However, because of plume meandering, steady 
wind shifts, lack of roadways and/or operator problems, it was not always possible 
to obtan the desired number and frequency of traverses. 

In the August 1979 field study, the remote sensors (i.e COSPEC and MARSS 
100) delineated the edges of the plume quite clearly. This was very helpful in 
determining the location to begin a new traverse. 

Also in the August 1979 study, the Mobile Air Monitoring Unit (MAMU) of 
the Monitoring and Instrumentation Development Unit carried out simultaneous 
stationary sampling along the traverse path of the near-field mobile vehicle. The 
purpose of this sampling was to assess the representativeness of the time-averaged 
cross-wind concentration profiles to stationary time-averaged concentrations. 
The MAMU was generally placed close to the center-line of the impingement zone 
of the plume and remained stationary for at least half an hour while the mobile 
unit traversed back and forth past it. Unfortunately, most measurement days were 
characterized by strong wind shifts and it was difficult to get long-term average 
data suitable for the desired analysis. 

Sampling in the MAMU was done for SO-, NO, N0 2 , NO x , H 2 S and O-j as 
well as surface level wind speed, wind direction, temperature, relative humidity 
and net solar radiation. Specifics of the instrumentation are available elsewhere 
(8). It should be noted that during the MAMU sampling program, the 
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plume often was subject to continuous wind shifts. Because of this, there was 
difficulty in obtaining enough long term average concentration data to do the 
desired comparison. 

Multipoint calibrations of the Sign - X instruments were carried out at least 
once every 2 days. The MAMU instrumentation was calibrated daily. 

During both field studies, a brief intercomparison test was carried out on 
the two Sign - X instruments. Both instruments, in their respective vehicles, were 
parked a few metres apart in a plume impingement area. The measured 
concentrations from the two units were compared and the results showed excellent 
agreement. 

2.2 METEOROLOGICAL MEASUREMENTS 

Meteorological measurements were carried out concurrently with the 
ground-level SO- concentration mesurements. Double theodolite pilot balloon- 
minisonde ascents were carried out for the measurements of wind direction, wind 
speed and temperature in the vertical. The ascents began early in the morning and 
continued until early evening. The interval between releases was roughly one hour 
in the 3une study and ^5 minutes in the August study. The continuous releases 
throughout the day were designed to provide information on the formation, growth 
and decay of the mixed layer as a function of time of day. 

The instrumentation consisted of two Warren-Knight theodolites, an Aero 
Aqua A A 3001 minisonde ground receiver and Aero Aqua balloons and minisonde 
transmitters. 
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The balloon release site was roughly 2 Km from the base of the INCO 381 m 
stack. The two theodolites were separated along baselines approximately 300 m 
long and aligned roughly perpendicular to the mean wind direction. The field 
procedures and data analysis methods were exactly the same as those used in the 
Sudbury Environmental Study Meteorological Field Program and are described in 
detail elswhere (9). 

Continuous measurements of total incoming solar radiation were also 
measured throughout both field studies. A Kipp and Zonen solarimeter of 
wavelength range 300 nm to 2.5 urn was used. It was permanently set up roughly 
2.5 km from the stack. The solarimeter output was used as a surrogate for surface 
heat flux in the model. 

2.3 STACK PARAMETER MEASUREMENTS 

INCO Ltd. operates continuous temperature, velocity and 50 2 
concentration monitors in the 381 m stack. The data are used in their Emission 
Reduction Program for the control of emissions under poor dispersion conditions 
(as required by a Ministry of the Environment Control Order). Hourly-averaged 
SO- emission rates were calculated by INCO as part of this system. Hourly- 
averaged values of buoyancy flux were calculated from the stack velocity and 
temperature measurments. 

All stack-related data presented in this report have been provided by INCO. 
Details of their measurement systems are available in other reports (10). 

2M HEAT FLUX MEASUREMENTS 

In the convective boundary layer, one of the scaling parameters controlling 
diffusion is the convective velocity scale, w # , defined in Equation 5 as 



w * = ( § H o V T,1/3 (m/s) 
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Here, g is the gravitational acceleration term, H is the surface kinematic heat 
flux, Z. is the average mixed layer height and T, the average mixed layer 
temperature. While the latter two variables are easily measured, the surface heat 
flux is not. Therefore, for this work the total incoming solar radiation, S, was used 
as a surrogate for heat flux according to: 



H = AS (W/m 2 ) (7) 

o 



where H is the maximum daily heat flux and A is a linear proportionality constant 
o 

typically varying from 0.1 to 0.55 depending on the surface. 



In developing the computer model, a value of A=0.3 was assumed. An 
experimental investigation of A was also carried out by making simultaneous 
measurements of heat flux and solar radiation during both field studies. A number 
of different methods were employed for these measurements- only one of which 
was successful. All are described here briefly. 

In June 1978, Dr. H. C. Martin of The Atmospheric Environment Service 
agreed to carry out micrometeorological measurements of heat flux and energy 
budget during the course of the field study. A micorometeorologial tower was 
erected in the valley north of Sudbury and measurements were taken on a half- 
hourly basis using eddy - correlation techniques. The experimental method, 
instrumentation and results of this work are described in detail elsewhere (12). 

Also in the June 1978 field study, attempts were made to indirectly measure 
surface heat flux using airborne measurments of temperature. The basis for this 
measurement is outlined in detail in Appendix 1. It can be summarized briefly by 
stating that the surface heat flux is related to the 
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standard deviation of temperature fluctuations, <ri measured in the convective 
boundary layer. Attempts were therefore made to measure c-£ from the same Piper 
Aztec aircaraft used in the concentration measurements program. The 
instrumentation for these measurements consisted of a Thermal Systems 
Incorporated fast-response hot-film anemometer system in the temperature 
measurement mode, a Hewlett-Packard Chart Recorder and a Hewlett-Packard 
3964A FM Recorder. The probe of the anemometer was mounted 20 cm beyond the 
nose of the aircraft. The control module and recorders were mounted inside the 
aircraft cabin. The TSI output was recorded both on strip chart and FM recorder - 
the former for visual observation and the latter for computer analysis. 

The measurements consisted of making two 30-50 Km cross-wind flights at 
a height of 100 m above ground level. These flights were made prior to the 
concentration profile measurments. 

Unfortunately this measurement program was unsuccessful. The 

temperature signal was continuously contaminated by RF interference and further 
use of the data was impossible. Spectral analysis of the recorded temperature 
signals indicated that many sources of interference were present and could not be 
eliminated. The method was abandoned mid-way through the study and no further 
reference will be made to the measurements. 

In the August 1979 study, another measurement method was used to infer 
the surface heat flux indirectly. Both the method and instrumentation were 
provided by Dr. B. Kerman of the Atmospheric Environment Service. Detailed 
descriptions of the theory, methodology and instrumentation are published 
elsewhere (13,1*0 but a short discussion is given here. 
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The instrumentation used was a laser scintillometer. It consisted of a 5 W 
He-Ne laser source and a remote photodetector. The post-detector signal 
processor produced a direct voltage output proportional to the variance of the 
logarithm of the received light intensity. This variance is related to the heat flux 
in a convective boundary layer and the algebraic relationship is known (see 
references). 

Early each morning the laser scintillometer was placed in a vacant field 
roughly 2.5 km from the 381 m stack. Approximately 50 m separated the laser 
from the detector. Proper alignment was obtained by using a theodolite mounting 
system. The data, in the form of output voltages, were recorded on an Esterline- 
Angus portable chart recorder. 

In analysing the data, a problem in the signal processing electronics was 
discovered. The effect of this problem on the measured data could not be 
quantified so further analysis was abandoned. 

3. DATA ANALYSIS AND RESULTS 

The meteorological and concentration data were analysed to produce results 
directly compatible with the computer model. This required that the 
meteorological data be analysed for mixed layer height, mixed layer wind speed, 
mixed layer temperature and incoming solar flux; the in-stack data for SC>2 
emission rate and buoyancy flux; and the concentration data for time-averaged 
centre-line concentrations. The individual methods of analysis and their results are 
discussed below. An analysis of the characteristics of the convective boundary 
layer in Sudbury is also presented. Note that all times referred to are Local 
Standard (Daylight Savings) Time. 
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3.1 CONCENTRATION DATA 

The cross-wind SO2 ground-level concentration profiles were analysed 
specifically to produce Eulerian average centre-line peak concentrations for direct 
comparison with model predictions. For the June 1978 study the data analysis was 
carried out by Moniteq Ltd. of Toronto. For the August 1979, study the analysis 
was done by a MOE summer student under the supervision of Ministry of the 
Environment staff as an undergraduate thesis for the University of Toronto. 

The method of analysis is a slightly modified version of a standard 
technique used in many other studies (5,6), It consisted of the following steps: 



(a) Each of the cross-wind SO^ profiles was screened to determine its usefulnes 
in the detailed analysis. Only blocks of data with a large number of 
consecutive traverses lasting greater than 20 minutes were accepted, 
although longer intervals were desired. Data were rejected if the profiles 
could not be completed (because the plume had passed beyond the end of the 
roadway) or if there was an insufficient number of traverses to provide a 
good average. 

(b) The chart records of the individual profiles were checked against 
topographical maps on which the traverse path had been marked. The 
ground reference points were identified and the known distances between 
reference points were checked against the chart record. 

(c) Individual SO2 profiles were grouped into subsets which were appropriate for 
time-averaging. These subsets required that the plume direction be 
relatively steady and the plume be continuously looping. Subset averaging 
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intervals ranged from a minimum of 19 minutes to a maximum of 135 
minutes for both studies. The minimum and maximum number of profiles 
averaged were 3 and 2<t respectively. The largest and longest set of 
traverses provided the most desirable and accurate data. 

(d) Each subset of profiles was analysed to produce Eulerian averages of gound 
level concentrations. An Eulerian average is defined as an average made 
with respect to a fixed reference frame (as opposed to a Lagrangian 
average which is made with respect to a moving reference frame). The 
analysis involved first digitizing the chart record profiles to produce 
concentration versus distance data along each traverse path. All of the 
individual profiles within a subset were then superimposed with respect to 
the ground reference points and the concentrations at each point along the 
traverse path were averaged. Figure 1(a) illustrates the individual profiles 
for one set of traverses and Figure 1(b) shows the resulting Eulerian 
average. Further details of the analysis method and the associated 
computer programs are given in References (15) and (16). Each subset 
Eulerian average was plotted and the magnitude and crosswind location of 
the peak concentration was obtained from the plot. This peak concentration 
was assumed to represent the centre-line concentration of the time- 
averaged plume and was directly comparable to the predicted 
concentrations from the model. 

(e) The distance and direction of the centre-line concentration from the INCO 
stack was determined from topographic maps with the traverse locations 
marked on them. 

(f) The Eulerian-average centre-line concentrations C(measured) and their 
respective distances from the stack were tabulated. These data appear in 
Columns 1 1 and 8 respectively in Tables 2 (a) and 2(b). . 
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3.2 METEOROLOGICAL DATA 

The pilot-balloon-minisonde data were analysed to produce values of mixed 
layer height, mixed layer wind speed, and mixed layer temperature as a function of 
time of day. 

To do this, the data were first analyzed to produce vertical profiles of wind 
speed, wind direction and temperature. The method of analysis and its associated 
computer program is described in another publication (9). Figure 2 shows a sample 
plot of one profile. The complete set for both field studies is presented in 
Appendix 2. 

From these profiles, a determination of the mixed layer height and mixed 
layer wind speed and direction was made. The method of determination was 
similar to that of Kaimal et al (1) and consisted of scanning the profile for the 
elevation at which a strong capping temperature inversion existed and/or the wind 
speed increased significantly after being relatively constant and/or the wind 
direction became constant after varying slightly. It was not necessary to use all 
of these criteria to determine the mixing height for every balloon release but, in 
most cases, more that one criterion was used. In some cases it was not possible to 
determine the mixed layer depth from the profiles because the measurements did 
not extend high enough. 

After the mixing height was determined, the representative mixed layer 
wind was calculated by resolving the North-South and East-West vector componets 
of the wind speed. These were calcualted by vertically integrating the component 
wind speeds to the top of the mixed layer. The resultant wind speed and direction 
therefore represented a vector-average, vertically-integrated, mixed layer wind. 
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Finally, the solarimeter data were digitized to produce hourly averages of 
solar flux. These are tabulated in Tables 3(a) and 3(b) . 

The daily variation of solar flux and mixing height is shown in Figures 3 to 
24. Days for which no concentration data were available are included in the 
figures if a reasonable amount of meteorological data was available. 



Once the solarimeter and mixed layer data were ready, further analyses 
were carried out to calculate the input data for the model. To do this, the solar 
flux, mixing height, and mixed layer wind speed were averaged over each interval 
for which Eulerian average concentration data were available. 

3.3 HEAT FLUX DATA 

The surface heat flux data obtained by Dr. H.C. Martin in June 1978 are 
presented in Reference (12). These data (obtained typically over 30 minute 
intervals) are also plotted in Figures 3 to 12 to show the diurnal variation of 
heat flux. Net radiation, Rn, data are also shown in these figures. 

3A SOURCE DATA 

The hourly-averaged in-stack temperature, velocity and SO- emission rate 
data were analysed for input to the model. This required estimating the actual 
time of emission of the portions of the plume sampled at ground level. This was 
done as follows: the average mixed layer wind speed and the distance from the 
stack for each set of traverses were determined. The time interval for each 
sampled portion of the plume to travel from the stack to its traverse location was 
calculated. This interval was then marked on a graph of emission rate, 
temperature and velocity versus time. An average value of each of these 
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parameters was determined from the graph for the estimated period of emission. 

For input to the model, the value of Q was taken directly from the stack 
data while the values of stack temperature and velocity were used to calculate the 
average buoyuancy flux from the stack, i.e. 

F = g V s R s 2 (T - Ta) (8) 

*s 



where F a buoyancy flux (m /K ) 

g = acceleration due to gravity (m/s ) 
V = stack gas velocity (m/s) 
R = radius of stack (m) 
T = stack gas temperature (K) 
Ta = ambient temperature at stack height (K) 



It should be noted that the INCO Emission Reduction Program often caused 
the emission rate of SO- to change quite significantly in a short period of time. As 
a result, the estimated emission rate averages presented here are subject to a 
great deal of variation and could seriously affect model validation. 

3.5 MODEL VALIDATION DATA 

The final data used to validate the fumigation model are presented in 
Tables 2(a) and 2(b). The date, start time, end time, averaging time and number of 
passes are shown for each averaging interval in Columns 1 to 5. Columns 6 and 7 
summarize the source data, i.e. the buoyancy flux, F and the SO- emission rate, 
Q. Column 8 gives the distance from the stack to the centre-line of the average 
profile. Column 9 lists the mixing height, Z. and Column 10 the 
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convective velocity scale, w # . Column 1 1 gives the mixed layer wind speed, u and 
the final column presents the measured, Eulerian average, centre-line SO ? 
concentrations. 

The accuracy of the data in Tables 2{a) and 2(b) is, in some cases, difficult 
to assess. The instantaneous individual stack parameters F and Q are accepted to 
be accurate within + 10% (verbal communication, INCO Ltd.). The distance from 
the stack should be accurate to within + 200 m and the height of the mixed layer 
should be within + 15% of the actual value (based on earlier comparisons of aircraft 
SO2 concentration and meteorological measurements). The estimates of 
convective scaling velocity have smaller uncertainties than any of the previous 
variables because of the one-third power dependence on the surface heat flux, the 
mixed layer height and the mixed layer temperature. Since it was not possible to 
measure surface heat flux and the adequacy of using solar flux as a surrogate is 
unknown, it was not possible to estimate the error in w # . The mixed layer wind 
speed had an estimated accuracy of + 20%. Finally, the Eulerian average 
concentration measurements are expected to be accurate to ± 25%. 

It should be pointed out that the model performance is directly dependent on 
the quality of the SO- emission rate data used in it. While the accuracy of a 
particular hourly average emission rate is + 10%, the uncertainties in the values 
shown in Tables 2(a) and 2(b) are much greater. This reflects the nature of the 
data and the requirements for analysing it. For example, on several occassions, the 
changes in SO_ emission rate with time were extremely rapid as the smelter 
responded to the requirements of INCO'S Emission Reduction Program. These 
changes were reflected in the hourly average SO_ emission rates measured from 
the in-stack monitors. However, the field measurements of ground level 
concentrations were always carried out over periods which did not correspond 
exactly to INCO'S hour-to-hour averaging period. It was 



- 19 - 

therefore necessary to estimate the emission rate for the measurement period by 
averaging the hourly - averaged data over the appropriate emission interval. 
Because of this, the impact of the rapidly changing emission rate on the measured 
ground level concentrations could be very significant. It is estimated that the 
uncertainty could be as high as + 40% in some cases. This can seriously affect 
model validation and should be carefully considered. It is apparent that a steady - 
state emission source is more desirable for studies of this nature. 

I*. DISCUSSION 

4.1 MODEL VALIDATION 

A comparison of the prediction of the short range, short term fumigation 
model with the field observations is presented elsewhere (4). The comparison was 
good and is discussed quantitatively in the reference. 

It should be noted that the centre-line concentration data used in the model 
validation have variable averaging times. As mentioned earlier, attempts were 
made to use the longest averaging time possible (within the limitations of the data) 
and the greatest number of traverses possible. To test the effects of the 
averaging time on the accuracy of the data and the model comparison, some of the 
August 1979 concentration data were analysed using different averaging intervals. 
That is, a given set of traverses was subdivided into several subsets which were 
averaged over their different time intervals. As expected, the comparison of 
model predictions with these concentration data showed a marked improvement 
with longer averaging times. A discussion of this is presented in detail in 
Reference (4). 

4.2 CHARACTERISTICS OF THE CONVECTIVE BOUNDARY LAYER 

The Dune and August intensive field studies provided an extensive data 
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base for determining the characteristics of the convective boundary layer in the 
Sudbury area. In general, these characteristics are very similar to those observed 
by Kaimal et al. (1) in Minnesota. 

In both the June and August field studies, the height of the convective 
boundary layer was generally well-defined as the height of a strong daytime 
capping inversion. Under this inversion the wind speed was relatively constant with 
height but accelerated quickly above it. The potential temperature within the CBL 
was also relatively constant. The wind direction typically varied with height 
within the mixed layer but became roughly constant above it. 

The diurnal cycle of the mixed layer followed a relatively consistent 

pattern. In the morning, the mixed layer developed rapidly as it responded to the 

increase in solar flux and the consequent increase in heat flux. This rapid growth 

generally occurred from 0900 to 1300 hours. At 1300 to 1400 hours, the mixed 

* 
layer reached its maximum height and remained relatively constant until 

approximately 1800 hours. 

Figure 25 illustrates this behaviour by showing the mixed layer height as a 
function of time for both field studies. The wide scatter in the data is a result of 
the changing metorological conditions which occurred during many of the days, 
such as rain events and overcast skies. However, the general pattern shows the 
rapid morning growth and the afternoon levelling off of the mixed layer. The 
maximum mixing height for both studies was 2300 m and the mid-day mixing height 
was typically around 1500 m. 

Figure 26 illustrates the diurnal variation of the total incoming solar flux. 
While a great deal of scatter exists within the data {again, due to varying 
meteorological conditions), a sinusoidal variation is apparent. The maximum daily 
solar flux was reached at 1300 to 1400 hours. After the maximum was 
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reached, the afternoon decline in solar flux had no apparent effect on the 
relatively constant mixed layer height. This was because the solar flux continued 
to generate sufficient convective energy to maintain the strong mixing in the CBL. 
Although there were no field data collected to investigate this in Sudbury, Kaimal 
et al. (1) suggest that the evening change in the mixing height is rather abrupt 
when compared to the morning growth. Their Minnesota data indicate that the flux 
and, therefore, the mixing height decreased very rapidly about an hour before 
sunset. The drop in heat flux propogated downward from the top of the CBL so 
that the surface heat flux was the last to turn negative. 

The late afternoon-early evening behaviour of the plumes from the INCO 
381 m stack and another 93 m stack (of Falconbridge Nickel Mines Ltd.) provided 
some evidence that the dissolution of the CBL in Sudbury also follows this pattern. 
For instance Kaimal et al. indicate that the heat flux turns negative from the top 
of the mixed layer down roughly to ground level at an hour before sunset. At this 
time in Sudbury, one can see the INCO 381 m stack plume begin looping with 
progessively smaller-scale loops while the Falconbridge plume continues to loop 
vigorously. This suggests that the vertical mixing is decreasing in intensity and the 
mixed layer height is continuously dropping in response to the decreasing heat flux. 
Within a very short time, the mixed layer height appears to have descended to 
just above the INCO 381 m stack. The plume then undergoes very small scale 
looping and most of it rises above the capping inversion over which it is trapped. 
However, the bottom portions of some loops are still caught in CBL and are 
brought to ground level by the remaining vertical downdrafts. Ground level 
impingements from these loops still occur but are less frequent, lower in 
concentration, and shorter in duration than afternoon ground level concentrations. 
While this is occurring, the Falconbridge plume continues looping. 
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Finally, near sunset, the CBL has descended to well below the effective 
height of the Superstack plume and no looping takes place. The plume at this time 
often takes on the appearance of a set of large discrete puffs with clear air in 
between. By sunset, both the INCO and Falconbridge plumes have stopped looping 
and convective activity has ended. 
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9 


11 


12 
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o 
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8-9 
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16 
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360 


49 
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311 


344 
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543 
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213 
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57 


410 


93 


41 


78 
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9-10 
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47 


492 


516 


311 


590 


57 


524 


549 


393 


541 
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756 


393 
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496 
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573 


129 


143 


398 


442 
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330 


51 
82 
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32 
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696 
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66 
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664 
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827 
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385 


688 
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20 


20 
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16 


20 
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FIGURE 1(a): INDIVIDUAL CROSS-WIND CONCENTRATION PROFILES 
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FIGURE Kb): EULERIAN AVERAGE CONCENTRATION PROFILE 
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Figure 2: Vertical Profile of Wind Speed, Direction and Potential Temperature 
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Figure 3: Diurnal Virlttlon of Meteorologies! Parameters 
(June 6, 1978) 
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Figure 4: Diurnal Variation of Meteorological Parameters 
(June 7. 1978) 
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Figure S: Diurnal Variation of Meteorological Parameters 
(June 8, 1978) 
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Figure 6: Diurnal Variation of Meteorological Parameters 
(June 9, 1978) 




Figure 7: Diurnal Variation of Meteorological Parameters 
(June 10, 1978) 
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Figure 8: Diurnal Variation of Meteorological Parameters 
(June 11, 1978) 
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Figure 9- Diurnal Variation oY'Mlieoroloo.lcal Parameters 
(June 13, 1978) 
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Figure 10: Diurnal Variation of Meteorological Parameters 
(June 14, 1978) 
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figure 11: Diurnal Variation of Meteorological Parameters 
(June IS, 197B) 
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Flaur* 12: Diurnal Variation of Meteorological Parameters 
(June lfi. 1978) 
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Figure 13: Diurnal Variation of Meteorological Parameters 
(June 20. 197B) 




Figure 14: Dlumal Variation of Meteorological Parameters 
(June 22, 1976) 
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F1gur* IS: Dfumil Variation of Meteorological Piramtttr* 
(June 23, 1578) 
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Figure 16: Diurnal Variation of Meteorological Parameters 
(June 24, 1976) 




Figure 17: Diurnal Variation of Meteorological Parameters 
(June 27, 1978) 
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Figure IS: Diurnal Variation of Meteorological Parameters 
(August 16, 1979) 



. R T -tQTtl MOM4M HIM 








- 
alMulf i,,IITI 


H*Bi*T*OM 










" 14 *HlttB UflM OIH* 










- 








- 


- 








- 


- 








- 


- 








- 


" 




'N* 




" 


" 


/k 




\ «i 


- 


■ 








- 


1 








■ 


/ 








- 


1 . , 1 1 J jf I 


J 1 L_ 


i i i 


J 1 1 L 


^TVj i i i 



MM 

MM 



i 

no* 

4*4 

tee 



■ova Q' t>** 



Figure 19: Diurnal Variation of Meteorological Parameters 
(August 17, 1979) 
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Figure ZO; Diurnal Variation of Meteorological Parameters 
(August 19, 1979) 
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Figure 21: Diurnal Variation of Meteorological Parameters 
(August 20, 1979) 
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Figure 22; Diurnal Variation of Meteorological Parameters 
(August 21, 1979) 
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Figure 23: Diurnal Variation of Meteorological Parameters 
(August 22. 1979) 



41 



. " T -10**1. (HCakMiM IOi.AH 






iteo 


- 


• 101 




IM« 




!•** 




HUM 




i»ft0 




ii« 




io<m 




»<J0 




• 00 




4C0 




*** 



Flour* 24: Diurnal Variation of Meteorological Ptruwters 
^ (August Z5, 1979) 
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Figure 25: Mixed Layer Height versus Time 

(June 1978 and August 1979 data) 
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Figure 26: Solir Flux versus Time Of Day 

{June 1978 and August 1979 data combined) 
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APPENDIX 1 

Using Temperature Measurements to 

Predict Surface Heat Flux 

A. Venkatram 

Under convective conditions the plume dispersion parameters ay and o z can De 
written as 

a y = a w t f y (t/T L ) ; a z = a w t f z (t/T L ) (1) 

In (1), a is the standard deviation vertical velocity fluctuation, t=x/u where u is 

the mean velocity, t is the time of plume travel, and T*l is the Lagrangian intergral 

time scale, At the short distances where high concentrations occur, it is possible 

to neglect the dependence of f v and f z on Tl # Thus, measurements of a and q 

j y z 

together with q and u will enable us to derive the functions fy and f z . Once a 

w y 

and <j are known, we can predict f if we can estimate . 
z y,z w 

Direct measurements of a at the levels at which elevated stacks emit pollutants 
poses serious experimental problems, and it is better to infer q through a more 
easily measured variable. Under convective conditions, this indirect measurement 
of (j is quite reliable because of the simple relationship between q and ot» the 
standard deviation of temperature fluctuations. We note that a- is a scalar 
variable which can be measured with relatively simple instrumentation. 

It can be shown that under unstable conditions a -r and q can be written as 
o w =A/g_HoZy/3; a T =B Ho 2/3/g_z\-l/3 < 2 ) 

where A = 1.33 and B = 1.35. In (2) H is the surface heat flux and T is the average 
temperature of the CBL. Form (2) we see that 



, w = i.i»/i t ^z\ 1 



11 (3) 
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It is worthwhile to point out that above Z >0.1Z,, o^ can be assumed to be 
constant (Zj is the mixed layer height). From (3) we see that measurement of a T 
with a fast response thermometer will enable use to estimate a^ and thus a and 

a z- 

We note from (2) that H is the only unknown in the formula for o w . It is 
clear than Oj allows us to compute H from the equation. 

H = (a T /B) 3/2 <gZ/T) 1/2 m 

Thus, in effect o"_ is a surrogate for H . In order to be able to predict o" w , we need 
to relate H computed from {*) to an operational available variable such as 
insolation (incoming solar flux). Simultaneous measurements of solar flux and Oj 
will then permit us to establish this relationship. Then, it is a simple matter to 
predict 0" and o and thus the pollutant concentration if we have a reasonable 
estimate of the incoming solar flux. 



It is useful to point out that a direct measurement of H with some sort of w • » 
device would be extremely useful in checking (4). However, direct measurement of 
H is less than completely reliable, and probably no more accurate than that 
inferred from (*). Thus, a comparison of H from (4) with H from w<0 ■ 
measurement would be more a test of consistency than accuracy of the proposed 
method. 
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Appendix 2 

June 1978 and August 1979 Meteorological Data 

Section A-2.1 : 

Plots of wind speed, wind direction and ambient temperature versus 

height are presented for June 1980. 

Section A-2.2: 

Plots of potential temperature versus height are presented consecutively 

for selected days in June 1978. 
Section A-2.3: 



Plots of wind speed, wind direction and potential temperature versus 
height are presented for August 1979. 
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